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We analyse the problem of a single mode field interacting with a pair of two level atoms. The 
atoms enter and exit the cavity at different times. Instead of using constant coupling, we use time 
dependent couplings which represent the spatial dependence of the mode. Although the system 
evolution is adiabatic for most of the time, a previously unstudied energy crossing plays a key 
role in the system dynamics when the atoms have a time delay. We show that conditional atom- 
cavity entanglement can be generated, while for large photon numbers the entangled system has 
a behaviour which can be mapped onto the single atom Jaynes-Cummings model. Exploring the 
I main features of this system we propose simple and fairly robust methods for entangling atoms 

, independently of the cavity, for quantum state mapping, and for implementing SWAP and C-NOT 

' gates with atomic qubits. 
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I. INTRODUCTION 



In recent years many authors have proposed different schemes, based on cavity QED systems, for entanghng atoms 
' and implementing quantum logic gates. Some of these proposals use single resonant interactions jlj or strongly detuned 
^p^: cavities [2;, ^ ^ for entangling atoms or carrying out logic gates. More elaborate schemes use decoherence-free spaces 
and continuous monitoring of the cavity decay for generating entangled states 0, S 0] ■ ^or these cases a no-photon 
^ , emission is associated with the generation of a maximally entangled state. In the optical cavity regime [s^, photon 
^ • polarisation measurements are used for entangling atoms [9f or implementing quantum logic gates [10, 11]. With these 
theoretical proposals the desired outcome is achieved by single-photon pulse scattering by an atom trapped inside an 
optical cavity jl2l | . 

As well as these theoretical proposals a number of experiments entangling Rydberg atoms in microcavities have 
[H [11, [S [H [II- EPR pairs [ij, [H, P ■ 



^ , been carried out [13|, [IJ, [15|, [16|, [17|, [18|. EPR pairs [IJ, [15|, [16| and entangled states of three quantum systems 
' [3, [13 were generated, a controlled phase gate and a C-NOT gate between a photon and an atom [3, [l3j were 
, implemented, a single photon non-demolition detection [l7| was perf ormed, and entanglement between an atom and 
■ a mesoscopic cavity field and measurement of decoherence effects [13] were carried out. Most of these experiments are 
[ based on the single atom Jaynes-Cummings ,19.] model where the atom and the cavity field, with zero or one photon, 
represent the qubits. The only exception is an experiment based on atomic collisions inside a strongly detuned cavity 
[3. [TgI for which the two atoms which collide inside the cavity are also the qubits. This experiment, in common with 
OO . the proposal we discuss here, is also relatively simple in that it does not require the use of Ramsey zones. So in 
' contrast to many of these proposals, in this paper we aim to obtain entangled atoms where the entanglement is not 
ultimately complicated by additional cavity entanglement. 

Among the theoretical proposals we also find the adiabatic, sequential passage of pairs of atoms through cavities 
^ , [13, [111. The scheme discussed by Marr et. al. [2C| exploits features similar to those appearing in STIRAP [2^. [23j. 

' The system adiabatically follows a dark state and a maximally entangled state between two atoms can be generated. 
" ■ In addition to this, Yong, Bruder and Sun [2l| proposed a scheme for generating entangled atomic states by sequential 
passage of atoms through a strongly detuned cavity. Using a time dependent Frohlich transformation, they derived 
an effective Hamiltonian and consider the preparation of entangled states with respect to the atomic velocities and 
the initial displacement between the atoms. We will see that the system we study here neither relies on the detection 
of lost photons, as Marr et. al. do [1^, nor requires accurate control of the delay between the atoms, as in the case 
of Yong et. al. [21]. 

We consider here as a key component of our scheme, the adiabatic limit for a system of two atoms resonantly coupled 
to a single mode field. Instead of the usual assumption of a constant coupling between the atoms and the field, we will 
utilise sequential time dependent couplings, which refiect the possibility that the two atoms enter the cavity mode at 
different times. From the analysis that follows, we show that the system is characterised by an energy crossing. This, 
and the symmetric structure of the adiabatic spectrum, make the system fairly robust; in an experiment one has to 
control a single mixing angle. Furthermore we show that a conditional atom-cavity entanglement can be generated and 
demonstrate a connection to the Jaynes-Cummings model which occurs when the cavity is highly excited. We also 
deduce the necessary conditions for using the adiabatic approximation (2,4^ , and discuss the feasibility of a potential 
experiment. Finally using our results we describe methods for entangling atoms and for mapping quantum states, and 
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we propose two setups for implementing a SWAP and a control-NOT gate (C-NOT). The proposed applications are 
fairly robust and simple to implement, whereas the operations could be relatively fast preventing potential failures 
due to decoherence effects. 

In contrast to previously proposed [E i, i, i, i, i, 0, i, [13, [H [ll or demonstrated methods [H, [11, [H, [H, [13, [Hi 

for entangling quantum systems or implementing quantum logic gates, here we consider a resonant two-atom, time- 
dependent Tavis-Cummings model, which differs from previous proposals either because they use a single atom 
interacting with a cavity, or because the interactions between the atom or atoms and the cavity are assumed to be 
constant in time. In addition to this, the main difference is that the resulting dynamics are different from what one 
expects and this is because of the existence of the energy crossing which is something new, and to our knowledge has 
never been reported before. 

The paper is organised as follows. In section [TT] we discuss our model and introduce the corresponding interaction 
Hamiltonian. After deriving the adiabatic states, we move to section Hill were we present and analyse our results. In 
section IIVI we discuss potential applications for quantum information processing and quantum state mapping. We 
conclude by summarising our results in section Ivl 



II. THE ATOM-CAVITY MODEL 



A. The Hamiltonian 



Our system consists of a pair of two-level atoms which resonantly interact with a single mode cavity field via time 
dependent couplings. In the interaction picture, and within the rotating wave approximation, the Hamiltonian reads 

2 

i=i 

The operators and a are the Bosonic creation-annihilation operators for the cavity mode, and <ji_ and cr^ are the 
lowering-raising operators for atom j. 

Both atoms are considered to traverse the cavity with the same speed v, following the same trajectory x{t), but 
atom 2 is delayed by 2At with respect to atom 1. The field spatial profile along the trajectory x{t) is assumed to be 
Gaussian 

Eix) = Eoexp{~^). (2) 
4a; Q 

The use of these assumptions results in Gaussian coupling pulses rij (r) : 

ry,(r)=ge~(-+^)^ r?,(r) = ge-(--^)^ (3) 
Here the dimensionless time r and the parameter S are defined in terms of the time width a = xq/v 

The width a could be seen as an average atom-cavity interaction time and later we will find that it can play the role 
of the interaction time in the Jaynes Cummings model . 

When deriving the Hamiltonian ([1]), we assumed that the atoms are moving sufficiently fast that the kinetic energy 
for each atom is greater than its coupling strength rjj{T). For this limit the atomic momentum operator p and the 
displacement operator x can be replaced by their classical counterparts mv and vt respectively. This requirement 
is easily satisfied with atomic beams, and we will also not expect the results of the variations to affect the form of 
the coupling shapes Tij^r). However, we not that in the opposite limit, when the atoms are slow, the atoms could be 
reflected without traversing the cavity psj . 

In what follows we consider a parameter regime in which the adiabatic theorem holds [23 | . Then if the Hamiltonian 
of the system is a smooth and slowly varying function of time, and the initial state of the system is one of the adiabatic 
states \'^j{ti)) of Hi{t), where 



Hj{t)\^,{t))^E,{t)\^,{t)), 



(5) 
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then at subsequent times t/, the system will be found in the adiabatic state \'^j{tf)) 

l*(i/)>=e-'*^(*^)|vl/,(i/)). (6) 

Eq. ([6]) is usually referred as the adiabatic approximation (23 | , where the dynamical phase acquired during the system 
evolution is 

<t>,it)^ J' dt'E.it'). (7) 

The conditions for applying this approximation will be discussed in detail later. We should also mention that the 
adiabatic theorem assumes that the Hamiltonian has a discrete spectrum, and for degeneracies the adiabatic approx- 
imation could fail. 

B. Adiabatic states 

The Hamiltonian Hj{t), Eq. ([1]), couples only bare states with the same number of total excitations. Thus our 
analysis is restricted to a subspace spanned by four states of the initial vector space, i.e. 

|n,eie2), |n + l,5ie2), \n + l,eigi), \n + 2,gig2). (8) 

The excited state of atom j is denoted \ej), \gj) is the ground state, and \n) refers to the field state with n excitations. 
The adiabatic energies are the roots of the characteristic polynomial P{E) 

P{E) ^ det\Ili{T) - E\, (9) 

where H/ is the matrix representation of the Hamiltonian in the subspace (jS]). From Eqs. ([T]) and ^ we have 

P{E) ^E^-E^{3 + 2n)(r/2 + -q^) + (1 + n){2 + n){ril - rjlf. (10) 

Because P{E) includes only even powers of E, the derivation of the corresponding roots is simple. The resulting 
adiabatic energies are 

Ei^2{t)=tE^{t), E3A^TE+iT), (11a) 



E^ir) . ^(3 + 2.)(,?(.)+,.^(.))±F„(.)^ ^^^^^ 

where the function Fnlj) is 

Fn{T) - ^ (77? (r) + ^l{T)f + 16(n + l)(n + 2WAtWAt). (12) 
With the energies in hand, Eq. ([5|) can be solved to give the adiabatic states, 

I^'i,2(t)) =A_(T)|n,eie2) +I?-(T)|n + 2,5i,g2) 

± (B_(r)|n + 1,5162) - C-{T)\n + 1,6132)), (13a) 
|^'3,4(t)) =A+{T)\n,eie2) +D+{T)\n + 2,gig2) 

± (B+(T)|n + 1,3162) -C+(T)|n + 1,6132))- (13b) 

The upper sign in Eqs. (jlSp is for the odd numbered states and the lower one is for the even ones. The coefficients 
■A±{t), B±{t),C±{t) and D±{t) are given in the appendix and Ref. [26^. For the purposes of our analysis it is 
sufficient to derive the limits of the adiabatic states for r — > ±cxd and r = 0. 
For the coupling functions ([3]), and for r = ±00, we assume that 

^lim = 0, lim f ^1 = 0. 
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FIG. 1: Top; The coupling functions ??i(r) (solid) and 7?2(t) (dashed) with respect to r. Bottom: The adiabatic energies -Ei(t) 
(solid), E2{t) (long dashed), -E3(r) (dashed) and E4,{t) (chain) with respect to r. The parameters are S = 1.0 and n = 0. 



Once these relations are used with Eqs. ()13a|) and (Il3bp . the adiabatic states for both limits are derived 

I f I", 6162) T 1?^-+ 1,5162) for T ^ -00 
^ 1 I", 6162) T 1?^-+ 1,6132) for T ^ 00 



I*i,2(t)) 



I*3,4(t)) - 



^ r |n + 2,3152) T |«+ 1,61.92) for 
>/2 I |„ + 2, .gi.92) T |n + 1, .9162) for 



T —00 

r — > 00 



(14a) 



(14b) 



From Eqs. (|14ap and (jl4bp we see that the adiabatic states, do not match individual bare states in both limits. 
This is well justified if we take into account the fact that the bare states are degenerate. Because of this even a small 
interaction is enough to lift the degeneracy and mix the bare states at the early stages of the system evolution. 

For T = 0, we have that Ei{0) — — E2{0). This can be seen either by putting rj^ = 772 in pT|) . or from figure [1] 
The corresponding degenerate states are also discontinuous functions of r 



|*i,2(t ^ 0=^)) = -(3n\n, 6162) + an\n + 2, 5132) ± 



sign(-r) , 
V2 ' 



The coefficients a„, /3„ and the state |— ) are 



1 



6 + 4n' 



6 + 4n' 



1 



|±) = -^|n + l)(|5ie2)±|eig2)). 



The remaining adiabatic states are continuous with respect to r 



1^-3,4(0)) = a„\n, 6162) + Pn\n + 2, gig2) T "^1 



(15) 

(16a) 
(16b) 

(17) 



Thus we see that the Hamiltonian Hj (r) has a temporal degeneracy. As already mentioned, when encountering such 
degeneracies the adiabatic approximation is expected to fail. If the degenerate states are coupled, then the system is 
likely to be found in a superposition of these states for r > 0, even if the coupling is relatively weak. A solution to 
this problem could be given by using a perturbation series [27| . Then, the resulting state vector will be expressed in 
terms of both adiabatic states, through the various terms appearing in the perturbation expansion. Here, instead of 
using a perturbation series, a simpler approach is considered based on an adiabatic elimination of two of the states. 
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III. SYSTEM DYNAMICS 
A. Energy crossing 

In order to predict the system evolution for t ~ we expand the Hamiltonian near the temporal degeneracy. From 
this expansion an effective Hamiltonian is obtained. The first step of this approach, is to change the current basis 
into that of the adiabatic states for r = 0^ , 

Hiir) - ^(-l)-'a;ir!-(T)|V,)(^,|+ ^(-l)^a;2f^+(r)|V,)(^,| 

J = l,2 i=3,4 

+%;^(l^2) - |^l))((V'3| + (V^4|) +h.c.. (18) 

The couplings Q± and the parameters a;i^2 are 

n±{T)=rj,iT)±Tj,{T), (19a) 



1 /(n + l)(n + 2) 1 



"^=2V 6 + 4n ' -2=2^6 + 4^' (19b) 
and l^jjj) — l^'j(O^). Then the Schrodinger equation in the new basis gives 

ici (r) = - (r) (^c^ici (r) + ^^^Zl±^^ , (20a) 
ic2(r)=0_(r)(^c.iC2(r) + ^^^^^±^), (20b) 
ic3(r) = - LU2n+{T)c3{T) + ^Mll(c2(r) - ci(t)), (20c) 

iC4(T) =^2f^+(T)c4(T) + ^^-ill (c2 (t) - Ci(r)). (20d) 

4W2 

Up to this point everything is exact. We now note that |f^+(T ^ 0)| 3> \^-{t ^ 0)|. This allows us to adiabatically 
eliminate {ips) and \ip4}, and derive an effective Hamiltonian. Setting c'a = = C4, and solving for 03^4 in terms of ci 
and C2, results in two differential equations 

iCl(T) =-Wir!_(T)ci(T), iC2(T) =Cc;ir!_(T)c2(T). (21) 

This equation corresponds to an effective Hamiltonian 

Heffir) = -4L0in-iT)a^ for r ~ 0, (22) 

where CTz is the Pauli matrix [2^ . 

The condition |f^+(T)| > |rj_(T)| holds for |t| < 0.25|(5|, with \d\ 1. For r < -0.25|5| we expect the adiabatic 
approximation to hold since the energy splittings between the adiabatic states are relatively large, see section FHI CI 
Assuming that the initial state of the system is |\['i(— oo)), then for r < — 0.25|(5| the system state reads 

l*(r)) = cxp (^-i dr'E,{r')^ |^Pi(r)) . (23) 

From Eq. ([22|) we see that |5'i(0~)) does not couple to \'^2{0 ))- Thus for t ~ the system state will be 

|*(0")) =exp ^-iy" dT'Ei{T')^\^i{0-)) T~0. (24) 
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This is true even for r = 0+, but since |^'i(0~)) = |'J'2(0"'")), Eq. ([15]), the state vector for r > will be 

l^-M) = exp (^-i J" dr' E^t') ~ i J\t' E2{t')^ |*2(t)) r > 0. (25) 

Thus for r = the system undergoes an energy crossing. The consequence of this is that if the system is initially 
prepared in cx))), it will end in state |4'2(oo)). Furthermore, since Ei{t) = —E2{—t) we can prove that 



/O noo 
dT'El{T')+ / dT'E2{T')=Q 
-oo JO 



(26) 



which means that the net dynamical phase for this pair of adiabatic states will be zero, i.e. 

|«'i(-^)) ^ 1*2(00)), |*2(-oo)) ^ |«'i(c»)). (27) 

In view of these results we can redefine the first two adiabatic states in order to properly incorporate the energy 
crossing so that 

Wiir)) - 0{-r) |*i(r)> + e{T) [^^(r)) (28a) 

l*2(r)) = e{~T) \^2{t)) + 9{r) \^,{t)) (28b) 



where 9{t) is the unit step function 



1 for r > 

e{T) = <( for T < (29) 
for T = 0. 
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The corresponding energies are 

E[(t) = -sign(T)ii;_(T), i?^(T) = sign(r)i?_(T). (30) 

To test the validity of Eqs. ((27l) . (|28ap and (j28bp we numerically solved the Schrodinger equation within a sym- 
metrical interval —Tq < t < Tq. In figure[2]the projection of Eq. ()28ap onto the state vector is plotted with respect to 
time. From this we see that the system adiabatically follows state \'^[{t)), with minor non-adiabatic effects. Despite 
this, the system ends in state |^'2(''"o)) in accordance with Eq. (j27p . In the lower graph the imaginary part of the 
projection (4'']^(ro)|5'(To)) is plotted with respect to the parameter ga. This is very small, in agreement with the 
expected zero value for the dynamical phase. The small deviations from zero can be related to minor non-adiabatic 
effects during the evolution. The choice made for the system parameters will be discussed later. 

B. Input-output in terms of the bare states 

For practical purposes, such as quantum information processing, one will encounter the bare states as input-output 
states instead of the adiabatic ones. Thus the results of section UlI Al must be expressed in terms of the bare states. 
Using Eqs. (I14ap and (|14bp the bare states are defined in terms of the adiabatic states for r —00. Then considering 
the limit for r ^ 00, taking into account P7)) and the dynamical phase acquired by the third and fourth adiabatic 
state, we have 

1?^, 6162) 1^,6162), (31a) 
|n-f 1,^162) ^ -|n -I- 1,61^2), (31b) 
\n + 1,61^2) -> cos{(t>n)\n + 1,5162) - isin((/)„)|n + 2,5152), (31c) 
\n + 2,5152) cos((/)„)|n -I- 2,5152) - ism{(f>n)\n + 1,5162), (31d) 



where the angle reads 



/oo 
drEiir). (32) 
-00 
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FIG. 2: Top: The projection |(5''i(r)|'I'(r))|^ with respect to time. The parameters are n — 0, S — 1.0 and ga = 30. The 
initial state was tq)), and |4'(t)) is the solution for the Schrodinger equation. Note the narrow range of the vertical axis. 

Bottom: The imaginary part of (^'i(ro)|4'(ro)) with respect to ga. The parameters are n = and S — 1.0. 



These four equations are enough to fully describe the system evolution. 

Eq. (j31b|) describes a complete energy transfer between the two atoms. This will happen without choosing special 
values for the system parameters provided we ensure the necessary conditions for adiabatic evolution. This robust 
energy transfer is a reminiscent of the STIRAP method [13, HI]. We also note that Eqs. picp and (j31dp describe a 
conditional entanglement of the second atom to the field mode. If atom 2 is not excited, then it will be entangled to 
the field mode. The exact form of the resulting entangled state is defined by 0„ for which we derive useful asymptotic 
expressions for various limits. 

The first two asymptotic expansions are those for (5 ^ 1 and S ^ 1. For the latter limit the mixing angle has a 
constant value proportional to \/n + 2: 



(f>n « 4:ga^/{n + 2)7r for 5 > 1. (33) 



This asymptotic expansion was derived by using an interpolation for the integrand in Eq. ()32|) . The interpolation 
hmction was the sum of two terms, the first is the limit of E^It) for r — cx) and the second one is the corresponding 
limit for t — > oo. For the other limit, 0„ is 



where 7„ is 



0„(i)»29ae-'Y2±i!*I fo, J«i (34) 



= Hiei±£±il. ,35) 



(3 + 2n)2 



This asymptotic expansion was derived with the Laplace method 2o|. The energy Ei{T) can be expressed in terms 
of a function w(t), such that £'4(t) — e~ For 6 the logarithm of 'w{t) can be substituted with a second 
order polynomial so that the integral results in Eq. (j34p . 

An interesting case of the asymptotic expression for is the one for n ^ 1. In this limit one can show that 
has the same dependence with respect to n as the mixing angle in the Jaynes-Cummings model [l^: 

4>n ~ 4:ga\/mr for n 3> 1. (36) 



This result suggests that for a large number of photons and with adiabatic evolution, we will have the same kind 
of dynamics as in the usual single atom Jaynes-Cummings model. Furthermore, this Jaynes-Cummings rotation is 
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FIG. 3: The matrix element Qlj^{T,5), Eq. (|37|l . with respect to r and S = 0.5, 1.0 and 1.25. The photon number was n — 
and the step dr — 0.01. 



conditional upon the state of the second atom. This could used for conditional operations in quantum information 
or for preparing field states with the use of conditional control. Of course in this limit the field dynamics will be the 
same as for the Jaynes-Cummings model [l9[ . 

The mixing angle 0„ can be tuned with velocity selection methods, so that a oc l/v. The sensitivity for these 
methods is of the order of A?; ~ 1% which corresponds to an error of the same order for the interaction time. 
The angle (/>„ can also be fine tuned with a Stark shift technique previously used in experiments with Rydberg atoms 
in microcavities [H [li, [H [H, [13, 111 • 

Finally, before discussing the necessary conditions for using the adiabatic approximation, we comment on the role 
of the first atom. The Eqs. (PTa|) - ((3Tdl) can been expressed in terms of a unitary operator U. This matrix can be 
factorised in terms of two unitary operators. The first one corresponds to a SWAP gate over the atomic basis with 
no effect on the field. The subsequent operation is a conditional rotation over the subspace of atom 2 and the field 
mode. The condition is that atom 1 must be in its ground state. Thus atom 1 could be seen as the "control qubit" 
during this sequence of unitary operations. 



C. Analysis of the adiabatic approximation 

Throughout sections IIII Al and IIII Bl we were assuming that the adiabatic approximation can be used. For this to 
be the case, the coupling between the adiabatic states must be smaller than the corresponding energy splitting |.24i] 

(M/,(r)|^|*,(r)) 

Qli ^) ^ N J « 5^' (37) 



\EI{t) - E'^irW 

where Hj{t) — Hj{T)/g and Ej{t) — Ej{T)/g. This choice of parametrisation makes Qn a function of n and S only, 
since the adiabatic states do not depend on g. Thus in order to quantify the conditions for adiabatic evolution we 
must be able to calculate the matrix elements of Q„. 

Using Eqs. (|13ap . (|13b|) and ([57)1 . and the fact that the coefficients A±{t), B±{t), C±(r) and D±{t) are real, it is 
easy to show that Q}^ and Q"^ are zero for any n and g and arbitrary time r. For the remaining QJf the derivation 
of analytic expressions is not simple. For this, we make use of numerical simulations, either by calculating or by 
solving the Schrodinger equation. 

The variation of Qn with respect to S and n, suggests that the adiabatic approximation can be used as long as (5 ~ 1 
and n is kept small. For n < 10 and ^ ^ 1, we have that Q„ ^ 1, figure [31 whereas for (5^1, Qn is substantially 
suppressed. Furthermore, we must ensure that ga 3> Qn at all times, i.e. 

gcr > max{(3„(T, (5)}. (38) 
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FIG. 4: The function (5e3(r) with respect to r and different ga: go = 5.0 (sohd), ga = 10.0 (long dashed) and ga — 20 
(dashed). The other parameters are n — and 5 = 1.0. 



For n < 10, the maximum of Q„ is of order 1 and Eq. ([55]) could be satisfied with ga > 10. If more excitations are to 
be added, then the lower bound is shifted towards higher values. Thus gcr must be increased in order to compensate 
for the increase in Qn- 

These conditions are confirmed by the results obtained after solving the Schrodinger equation. Using as input the 
adiabatic states, we integrate the Schrodinger equation for a time interval |t| < Tq. Then to quantify the non-adiabatic 
effects, and understand how they can be suppressed by proper tuning of the physical parameters, we use the function 

Se,iT)^\l-\{^,{r)\^iT))f\. (39) 

What we find from these simulations is that the delay parameter 5 must be of the order of unity, i.e. 1 < S < 1.25. 
Furthermore the coupling amplitude g and the interaction time a, must be large enough so that their product exceeds 
a lower bound which is a function of the photon number n. For n = 0, this lower bound is approximately equal to 10, 
figure m and ga 3> 10. For larger n, the bound increases, figure O and the coupling amplitude and the interaction 
time must become larger in order for the bound to be satisfied. 

The violation of the adiabatic approximation with increasing photon number is also shown in figure[6l For increasing 
values of n, and assuming that the product ga is constant, we see that the adiabatic approximation will fail. Thus 
this approximation has a dependence with respect to n. Then if we want the system to evolve adiabatically for larger 
photon numbers, a stronger coupling must be chosen. 

Physically the requirement ^ ~ 1, means that while the atoms must both be in the cavity at some point, they must 
enter and exit the cavity at different times, figure [T] This demands control over the delay time and the velocities of 
the atoms. Furthermore, condition (|38p requires strong couplings g, and suppression of decoherence to ensure longer 
interaction times a. 

A typical optical cavity [s^l has a lifetime of the order of 40ns while the coupling strength is approximately 200MHz. 
This means that the dimensionless product ga cannot exceed 10. For a micromaser cavity [sj, with quality factor of 
the order of 10* the photon lifetime is approximately 160/is, with the effective interaction time being equal to 20^s. 
The coupling is about 150kHz which gives ga ~ 3. For even higher quality factors, e.g. lO^*', the photon lifetime is 
of the order of 0.1s which enables interaction times of the order of ct ~ 100/is [32]. Then with the same coupling, 
g ~ 150kHz, we have that ga ^ 15 which could be enough for observing adiabatic evolution with ti ^ 1. This results 
in sufficiently fast atoms to overcome the barrier due to the coupling, e.g. for Rb atoms v ~ lOOm/s Thus a 
high quality micromaser cavity appears to be a good candidate for realising the adiabatic evolution for low photon 
numbers. Moving towards higher n would require either higher quality factors or stronger fields. 
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FIG. 5: The function 5e3{T) with respect to r and different n: n = (solid), n = 5 (long dashed) and n — 10 (dashed). The 
parameters are ga = 30 and S — 1.0. 



IV. APPLICATIONS 



A. Atomic entanglement 



We have already seen in section UlI Bl that a single passage of the atoms through the cavity can entangle the atoms 
to the cavity. For atomic entanglement, we would like the cavity to be disentangled. In order to do this we must 
prepare the atoms in superpositions of the form aj\gj) + j3je^^^ \ej) while the cavity is kept empty 

|*m> = |0) ® (ailgi) +/3ie*^i|ei)) ® (a2|52> + /32e*'^ jea)) . (40) 

The coefficients aj and (3j are real and aj + (3j — 1. After the passage of the atoms through the cavity and assuming, 
for now, that 0_i is arbitrary, the output state will be 

|*out> = |0) ® l^-en) - i/?ia2 sin((/._i)e"'i 1 1) ® I5152), (41) 

where the atomic state |^cn) reads 

|*cn) = ai\g2) {a2\9i) - I32e'''\ei)) 

+/?ie''i |e2) (a2 cos(0_i)|gi) + f32e''^ |ei)) . (42) 

Thus, if the mixing angle 0_i is not equal to (2m + l)7r, where m = 0, 1, 2 . . . , the atomic state l^'cn) will be an 
entangled state of the two atoms. The probability Pen for this state is 

Pen = l-/3?a2sin2(0„i)- (43) 

If we want to generate a maximally entangled state, then we must ensure that Pen = 1 by choosing 0_i = 27r. 
Furthermore aj and /3j must be l/\/2 and the phase angle 62 is zero. Then the resulting state is 

l*c„) = I (I52) (|<?i> - \ei)) + e'"^ |e2) (I51) + |ei))) , (44) 
and a Hadamard rotation H over the subspace of atom 1 33"| must be applied 

i/|5i) = ^(|gi> + |ei)), i/|ei) = i=(|5i)-|ei)), (45) 



to get the maximally entangled state 



Ivfen) - ^ (e'^M3ie2) + |eig2)) . (46) 
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FIG. 6: The function 5e3{T) for r ^ oo with respect to n. Top: ga — 10 (solid), ga — 20 (long dashed) and ga = 30 (dashed). 
Bottom: ga — 50. For both plots we have S = 1.2. 



B. State mapping 



Assume now that the two atoms are prepared in the following superposition 

|gi)(a|g2)+/3|e2)), (47) 

and the cavity is in the vacuum state, then the state of the system after the passage of the atoms through the cavity 
will be 

|0,52)(a|5i)-/3|ei)). (48) 

Applying a single rotation on atom 1, |ei) — > — |ei), the final state of atom 1 reads 

a|3i)+/3|ei). (49) 

Thus the state of atom 2 can be mapped into atom 1. The whole scheme is simple and robust since there is no need 
for exact control over the interaction time, in contrast to other proposals [l[. 

In addition to this scheme, one can implement state mapping between the atoms and the cavity, and vice versa. 
If the cavity is prepared in the state a|0) + /3|1), the two atoms are in their ground states, and the mixing angle is 
— tt/2, then the resulting state for the second atom will be a\g2) — i/3|e2) whereas the cavity and the first atom 
are not excited. If the single qubit rotation |e2) i|e2) is applied, then the resulting state for atom 2 will be 

a|32)+/3|e2). (50) 

If now the system is prepared in the state |0,g2) + /3|ei)), and we apply the transformation |ei) —^ i|ei), and 

send the atoms through the cavity with = 7r/2, then the state of atom 1 will be mapped onto the cavity mode, 
i.e. the system state becomes 

(a|0)+/3|l))|5i,<?2). (51) 



C. SWAP and C-NOT gates 



We conclude the analysis of the current section by presenting two setups for implementing a SWAP and a C-NOT 
gate. The first can be realised using an empty cavity and choosing 0_i = tt. From Eqs. (|31ap - (|31cp we have 

|eie2) \e1e2), \g1e2) -\e1g2), 
\91g2) \91g2), \e1g2) -15162), 
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with the cavity remaining empty. Thus the model in hand could be used to implement a SWAP gate, without using 
any additional components. 

In contrast to this, the realisation of a C-NOT gate requires the use of additional components. We consider a row 
of two empty cavities, with the angle 4>-i being 2tt and tt respectively. If a pair of atoms crosses through these two 
cavities, then the output is 

\e1e2) \e1e2), \g1e2) \g1e2), (52a) 
15152)^15152), |ei52) ^ -|ei52)- (52b) 

The use of an array of two cavities is allowed since neither cavity is excited after the atoms have passed through. In 
addition to this, attention must be paid, so that both atoms exit the first cavity, before entering the second one. 

Eqs. (|52ap and (j52b[) correspond to a phase gate. This can be combined with two Hadamard gates to produce a 
C-NOT gate. For the current setup, the Hadamard rotation is performed over the qubit space represented by atom 1, 
Eq. (|45|) . and is applied before and after the atoms have crossed the cavities. Combining Eqs. ([45]) . (I52a|) and (|52b|) 
gives the following input-output table 

\e1e2) leiCa), 
15162) \91e2), 
\e1g2) I5152), 
I5152) \e1g2), 

which corresponds to a C-NOT gate, with atom 2 being the control qubit and atom 1 the target qubit. 

With the exception of the state mapping between atoms, which is a robust scheme, the other proposed applications 
in this section could be characterised as fairly robust. The fidelities for these applications can exceed 99% with proper 
control of the physical parameters such as interaction time a and delay time At. For example, in figure [7] the fidelity 
for a maximally entangled state is plotted with respect to Act for various errors in the delay time. From this we see 
that the fidelity can be greater than 97% if the error in a is below 5% while the corresponding error for the delay 
time can be as high as 10%. With better tuning of the interaction time and after suppressing the error Act below 2% 
the fidelity can exceed 99%. 

In general for all the given applications the corresponding fidelity is less sensitive to errors in the delay time in 
contrast to errors in the interaction time a. This means that the tuning of the interaction time has to be more accurate 
than the corresponding tuning for the delay time. Furthermore, for some of the applications, such as mapping the 
cavity state (|0) -I- |l))/\/2 onto atom 2, the threshold for the error in a could be as high as 10% with a resulting 
fidehty equal to 99%. 

Decoherence effects related to spontaneous emission from the atoms and the dissipation of photons from the cavity, 
are detrimental to the proposed applications. The adiabatic condition, Eq. (|38p . suggests that the mean interaction 
time can be substantially decreased while increasing the coupling strength g to ensure adiabatic evolution. This 
will allow fast implementation of the proposed applications, reducing the probability for potential failures due to 
decoherence. 

Finally, we note that the implementation of the single atom rotations can be carried out with classical EM fields 
[S^ l . For the rotations to be successful, we must ensure that the spatial displacement between the atoms is sufficiently 
large. This will prevent interactions between the fields and both atoms when only one atom must couple to the EM 
fields. 



V. CONCLUSIONS 

In this paper, we consider the situation where two atoms interact with a single mode cavity field via sequential 
time dependent couplings. We focused on the adiabatic limit, but were also able to show that an energy crossing 
takes place in the vicinity of a temporal degeneracy. This takes place because the degenerate adiabatic states do not 
couple to each other near the degeneracy point. This effect, together with the symmetry properties of the adiabatic 
spectrum, makes the system fairly robust. The only parameter which must be controlled is a mixing angle This 
angle defines the degree of conditional entanglement between the atoms and the cavity. In the limit of large photon 
number, the mixing angle has the same dependence with respect to the photon number as in the usual single atom 
Jaynes-Cummings model. Thus in this limit we could say that the system behaves like a conditional Jaynes-Cummings 
system. 

A significant difference frorn previously proposed theoretical schemes and experiments carried out [J, [2j S 0, S S 
0, H, [13, [m [13, [H, [13, [IB [3 [S [3 J is the energy crossing which is seen in figure [TJ This feature enables the robust 
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FIG. 7: The fidelity for a maximally entangled state with respect to variations in the interaction time Act and different variations 
in the delay time: (solid), 0.05At (dashed) and O.lAt (dot). The parameters a and At correspond to S = 1.0, t/i-i — 2mn 
and 01 = 0. 

disentanglement of the atoms from the cavity, after the evolution of the system is completed. Based on this, we have 
proposed methods for entangling atoms and mapping quantum states between atoms and cavities and experimental 
setups for implementing a SWAP and a C-NOT gate. The proposed schemes are fairly robust, with fidelities up to 
99%. This robustness is especially evidenced in the delay time since errors of up to 10% have minor effect on the 
fidelities of the proposed applications. Furthermore, the fact that the interaction time can be substantially decreased 
makes these applications less sensitive to decoherence effects. 

What makes this system fairly robust is the symmetric structure of the adiabatic spectrum and the energy cross- 
ing. Both of these features are sensitive to variations of the detuning, and to changes of the ratio of the coupling 
amplitudes. If the coupling amplitudes have a ratio that is slightly different from unity, then the system evolution 
will be significantly different. This should be expected since a second mixing angle appears in Eqs. pip , and one 
must be able to control the ratio gi/52 in order to have the desired output. Furthermore, the existence of a detuning 
between the atoms and the field will lift the temporal degeneracy, preventing the energy crossing from occurring. The 
detuning also affects the structure of the spectrum and its symmetry properties. 

The system is vulnerable to decoherence effects. Since we are interested in the adiabatic limit, the interaction time 
is expected to be large. From our analysis we conclude that the adiabatic condition allows shorter interaction times 
by simply increasing the coupling strength. In this way a faster process can be realised, reducing potential failures 
due to decoherence. Based on current experiments, either with optical or micromaser cavities, we find that a high 
quality micromaser cavity could be used for realising the proposed applications. 
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APPENDIX A: ADIABATIC STATES 



The adiabatic states are derived after solving Eq. ([5]) . The normahsed adiabatic states are given in p3a|) and (pb]) 
where the coefficients A±{t), B±{t), C±{t) and D±{t) are 



A + 4(n+l)(n + 2)7727^2 



S±(t) 



2r,,V^rT2 (i?|(T) + {n + - r,l)) 



(Alb) 



C±{t) 



jJ±(T)^±(T)r/i(3 + 2n) 
V^rT2(i?^(r) + (n + l)(r,2-r;2)) 



(Ale) 



Similar expressions were previously derived in the context of pair effects in single atom micromasers [26j . 



[1] S. B. Zheng, Phys. Rev. A 71, 062335 (2005). 

[2] S. B. Zheng and G.-C. Guo, Phys. Rev. Lett. 85, 2392 (2000). 

[3] E. Jane, M. B. Plenio, and D. Jonathan, Phys. Rev. A 65, 050302 (2002). 

[4] L. You, X. X. Yi, and X. H. Su, Phys. Rev. A 67, 032308 (2003). 

[5] M. B. Plenio, S. F. Huelga, A. Beige, and P. L. Knight, Phys. Rev. A 59, 2468 (1999). 

[6] A. Beige, S. Bose, D. Braun, S. Huelga, P. Knight, M. Plenio, and V. Vedral, Journal of Modern Optics 47, 2583 (20 
November 2000). 

[7] A. Beige, D. Braun, B. Tregenna, and P. L. Knight, Phys. Rev. Lett. 85, 1762 (2000). 

[8] R. Miller, T. E. Northup, K. M. Birnbaum, A. D. B. A Boca, and H. J. Kimble, Journal of Physics B: Atomic, Molecular 
and Optical Physics 38, S551 (2005). 

[9] L. M. Duan and H. J. Kimble, Phys. Rev. Lett. 90, 253601 (2003). 
[10] L. M. Duan and H. J. Kimble, Phys. Rev. Lett 92, 127902 (2004). 
[11] L. M. Duan, B. Wang, and H. J. Kimble, Phys. Rev. A 72, 032333 (2005). 
[12] L. M. Duan, A. Kuzmich, and H. J. Kimble, Phys. Rev. A 67, 032305 (2003). 
[13] J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. Phys. 73, 565 (2001). 

[14] G. Nogues, A. Rauschenbeutel, S. Osnaghi, M. Brune, J. M. Raimond, and S. Haroche, Nature 400, 239 (1999). 
[15] E. Hagley, X. Maitre, G. Nogues, C. Wunderlich, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 79, 1 (1997). 
[16] S. Osnaghi, P. Bertet, A. Auffeves, P. Maioh, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 87, 037902 
(2001). 

[17] A. Rauschenbeutel, G. Nogues, S. Osnaghi, P. Bertet, M. Brune, J.-M. Raimond, and S. Haroche, Science 288, 2024 
(2000). 

[18] A. Rauschenbeutel, G. Nogues, S. Osnaghi, P. Bertet, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 83, 
5166 (1999). 

[19] M. O. Scully and M. S. Zubairy, Quantum Optics (Cambridge University Press, Cambridge, 2002). 
[20] C. Marr, A. Beige, and G. Rempe, Phys. Rev. A 68, 033817 (2003). 
[21] L. Yong, C. Bruder, and C. P. Sun, Phys. Rev. A 75, 032302 (2007). 

[22] K. Bergmann and B. W. Shore, in Molecular Dynamics and Spectroscopy by Stimulated Emission Pumping, edited by H. L. 

Dai and R. W. Field (World Scientific, Singapore, 1995), chap. 9, pp. 315-73. 
[23] K. Bergmann, H. Theuer, and B. W. Shore, Rev. Mod. Phys. 70, 1003 (1998). 
[24] A. Messiah, Quantum Mechanics (Dover Publications, New York, 1999). 
[25] G. M. Meyer, M. O. Scully, and H. Wahher, Phys. Rev. A 56, 4142 (1997). 
[26] S. Mahmood and M. S. Zubairy, Phys. Rev. A 35, 425 (1987). 
[27] G. A. Hagedorn, Annals of Physics 196, 278 (1989). 

[28] R. Wong, Asymptotic approximations of integrals (Society for Industrial and Applied Mathematics, Philadelphia, 2001). 
[29] B. T. H. Varcoe, private communication. 

[30] A. Boca, R. Miller, K. M. Birnbaum, A. D. Boozer, J. McKeever, and H. J. Kimble, Phys. Rev. Lett. 93, 0233603 (2004). 



15 



[31] M. Bruno, E. Hagley, J. Dreyer, X. Maitre, A. Maali, C. Wunderlich, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 

77, 4887 (1996). 

[32] B. T. H. Varcoe, S. Brattkc, and H. Walthcr, New Journal of Physics 6, 97 (2004). 

[33] M. A. Nielsen and 1. L. Chuang, Quantum computation and quantum information (Cambridge University Press, Cambridge, 
2000). 

[34] F. Yamaguchi, P. Milman, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. A 66, 010302 (2002). 



